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ABSTRACT 

The central star-forming regions in three blue compact dwarf galaxies (He 2-10, NGC 5253, 
and II Zw 40) were observed in the 340 GHz (880 [xm) band at ^ 5 arcsec resolution with the 
Submillimetre Array (SMA). Continuum emission associated with the central star-forming 
complex was detected in all these galaxies. The SMA 880 |J.m flux is decomposed into free- 
free emission and dust emission by using centimetre-wavelength data in the literature. We 
find that free-free emission contributes half or more of the SMA 880 [im flux in the central 
starbursts in those three galaxies. In spite of the dominance of free-free emission at 880 i.im, 
the radio-to-far infrared (FIR) ratios in the central star-forming regions are not significantly 
higher than those of the entire systems, showing the robustness of radio-FIR relation. Based 
on the robustness of the radio-FIR relation, we argue that the free-free fraction in the 880 \xm 
emission is regulated by the dust temperature. We also analyze the CO (J ~ 3-2) emission 
data. We find that CO is a good tracer of the total gas mass in solar-metallicity object He 
2-10. Low-metallicity objects, NGC 5253 and II Zw 40, have apparently high star formation 
efficiencies; however, this may be an artifact of significant dissociation of CO in the low- 
metallicity environments. We also point out a potential underestimate of dust mass, since the 
dust traced by emission is biased to the most luminous high-temperature regions, particularly 
when a system hosts a compact star-forming region where the dust temperature is high. 

Key words: dust, extinction — galaxies: dwarf — galaxies: evolution — galaxies: individual 
(He 2-10, NGC 5253, II Zw 40) — H II regions — submilHmetre: galaxies 



1 INTRODUCTION 



The early stage of galaxy evolution can be characterized by 
the following two properties: a poor metal abundance and a 
, rich gas content. Although galaxies at the early evolutionary 
stages exist at high redshift, nearby blue compact dwarf galax- 
ies (BCDs) are a unique category of galaxies that have those 
two properties in the nearb y Un iverse (Sargent & Searle 1970; 

Ivan Zee. Skillman. & Salzerl il998l : iKunth & Ostlirj |2000|) . Some 
BCDs are also experiencing the most active cla ss of star formation 
with the formation of super st ar clusters (SSCs) (iTumer et al.ll99^ : 

iKobulnickv & JohnsorJ ll999'). This kind of intense star formation 
may also provide a relevant laboratory to understand the condi- 
tions in which high-redshift galaxies experience their first starburst 
episode. 

Such intense star formation as seen in BCDs is occurring in 
dense and compact regions and thus can only be traced with op- 
tically thin star for mation indicators, such as far-infrared (FIR) 
dust l uminosity (e.g. lKennicutl|[l993 : llnoue, Hirashita. & Kamaval 
l2000h and radio luminosity (thermal free-free emission from H II 
regions plus non-thermal synchrotron emission from accelerated 
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electrons) dCondonll 19921) . These two luminosi ties are strongly cor- 
related in nearby star-forming galaxies (e.g. Ide Jong et al.lll98^ : 
iHelou, Soifer. & Rowan-Robinson|[l98^. 

In our previous paper, iHirashital fcOllL hereafter Hll), we 
observed a w ell studied BCD, I I Zw 40, with the Submillimetre 
Array (SMA; IHo, Moran. & Lol 12004) . In this paper, we add two 
BCDs, He 2-10 and NGC 5253, both of which also host high star 
formation activities likely to be associated with the formation of 
SSCs in the centre ( lGoriian|[T995 Ijohnson et al. 2000). The radio 
emission at wavelengths > a few cm is optically thick for free- 
free absorption, supporting the existenc e of young com pact regions 
before the expansion of H I I regions (Turner. Ho. & Beck 199^ 
IKobulnickv & Johnsonl 1 19991 : Ijohnson & Kobulnick^i 12003 )EI All 
these BCDs are also classified as Wolf-Rayet galaxies: the Wolf- 
Rayet feat ure indicates that the typical age of the current starburst i s 
a few Myr l iVacca & Contill992[ : lL6Dez-Sanchez & Estebarj20"lol) . 

Some of the conclusions derived from II Zw 40 by Hll can 
be generalized by increasing the sample. Hll argues that free- 
free dominated submillimetre emission can be a characteristics of 

^ In this paper, we focus on the wavelength range 2 cm), where the 
radiation is optically thin for free-free absorption as far as our sample is 



concerned (Sections l2. 1H2.3 
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young active star formation. Hll has clarified that, if we focus on 
the central star-forming region in II Zw 40, free-free emission dom- 
inates even at 880 [ivci. Consi dering that globafl su bmm luminosity 
is usually dominated by dust jGalliano et al ]|2005l) . the dominance 
of free-free emission at 880 ^m is a special characteristics, which 
should be examined with a larger sample. 

By increasing the sample, we can also explore the metal- 
licity dependence of various properties. The metallicity of 
the interstellar medium is a fundamental quantity character- 
izing the galaxy evolution because it reflects the enrichment 
in heavy elements by stellar generations. Some metallicity ef- 
fects are expected in the formation of dense star-forming re- 
gions: (i) Dust-to-gas ratio in low metallicity objects is gen- 
erally low /Schmidt &. Boiler 'l993|; iLisenfeld & Ferraral 1 19981: 
iHirashita , Taiiri, & Kamava 2002), which implies that the star- 
forming regions are less embedded by dust in low-metallicity ob- 
jects than in high-metallicity objects; and (ii) in a low-metallicity 
(i.e. dust-poor) environment, star formation is expected to be less 
efficient because of less shielding of ult raviolet (UV ) heating pho- 
tons by dust and molecular hydrogen ( Hirashita & Ferrara 2002; 
lYamasawa e"tar]|20Ill; [Onedin & Kravtsovll201 ill Thus. we will 
address the metallicity dependence of dust abundance and star for- 
mation efficiency by investigating gas and dust emission in BCDs. 

In observing intense starbursts associated with the formation 
of SSCs, a high spatial resolution is crucial to spot the young star- 
forming component. While high-resolution data of the central star- 
forming regions in BCDs are available through various radio inter- 
ferometric observations with arcsecond resolutions (1 arc sec cor- 
responds to 51 pc for He 2-10 and II Zw 40, and 18 pc for NGC 
5253), information of FIR-submm dust emission on such a small 
scale was lacking for BCDs. We thus performed SMA observa- 
tions of a few BCDs to resolve their submm emission. Shorter 
wavelengths such as mid-infrared can achieve almost arcsecond 
resolutions by single-dish telescope facilit ies such as the Sp itzer 
Space Telescope Infrared Array Camera ( Fazio et all |2004|) and 



AKARl Infrared Camera dOnaka et al.ll2007h . but the mid-infrared 
emission is dom inated by stochastically heated very small grains 
( iDraine & And erson 1985), which are not representative of the to- 
tal dust amount (although it is empirically known that the mid- 
infrared l uminosity is we ll correla ted with the total dust luminosity 
in BCDs IWu et alj|2008l ; see also lTakeuchi et alj|2005h ). 

This paper is organized as follows. We explain the observa- 
tions and the data analysis in Section |2] In Section [3] we derive 
basic quantities related to star formation and submm emission, and 
discuss radio-FIR relation. We also analyze CO(3-2) emission and 
derive the molecular gas mass, which is used to estimate the dust- 
to-gas ratio and the star formation efficiency. After discussing our 
observational results in Section |4] we conclude in Section |5] 



2 OBSERVATIONS AND DATA ANALYSIS 

We selected nearby BCDs which host young compact starbursts. 
Searching for objects feasible for the sky coverage and sensi- 
tivity of SMA, we chose He 2-10 and NGC 5253. They have 
similar properties to our previous object, II Zw 40, in Hll: (i) 
The ages of the central starbursts are y oung (< a few Myr) 
as indicated by the Wolf-Rayet feature jVacca & Contil 1 1 9921 ; 



'Lopez-Sanchez & Esteban"201(f) and the stellar spectral synthe- 
sis models (Chandar et al. 2005); and (ii) they host optically thick 
compact free-free emission at wavelengths > a few cm, indicating 
an intense starburst in a dense compact region ( Turner et"ani 19981 ; 
iKobulnickv & JohnsoiJll999l ; |johnson & Kobulnick'sll2003l ). Since 
He 2-10 has a higher metallicity than the other two galaxies, we 
may be able to obtain a hint for metallicity effects (Section [3] see 
Table[T]for the observed metallicity values). The distances (D) and 
velocities (relative to the local standard of rest; Vlsr) adopted in 
this paper are listed in Table [T| 

The SMA observations of He 2-10 and NGC 5253 were car- 
ried out in the 340 GHz (880 \im) band on 2012 Januaiy 11 in 
the subcompact configuration. Seven antennas were used with pro- 
jected antenna separations between 9.5 and 45 m. The receivers 
have two sidebands, the lower and upper sidebands, which cov- 
ered the frequency ranges from 330.8 to 334.8 GHz, and from 
342.8 to 346.9 GHz, respectively. The visibility data were cali- 
brated with the MIR package. As a flux calibrator we used Cal- 
listo (with an adopted flux of 14.5 Jy) for He 2-10, and Titan 
(with an adopted flux of 2.36 Jy) for NGC 5253. We used quasars 
J0730— 1 16 and J1337— 129 as amplitude and phase calibrators for 
He 2-10 and NGC 5253, respectively. We adopted quasar 3C279 
as a band pass calibrator for both objects. In generating the contin- 
uum data, we excluded the chunk containing the CO(3-2) emission 
at rest 345.796 GHz. The calibrated visibility data were imaged and 
CLEANed with the MIRIAD package. The synthesized beam has a 
full width at half-maximum (FWHM) of 5.7 arcsec x 4.4 arcsec 
(290 pc X 220 pc) with a major axis position angle of 41° for He 
2-10 and a FWHM of 6.7 arcsec x 4.1 arcsec (120 pc x 74 pc) 
with a major axis position angle of 33° for NGC 5253. The largest 
angular scale sampled by this observation is 19 arcsec. We also use 
the SMA subcompact data of II Zw 40 in Hll (see HI I for details). 

Figure [T] shows the obtained continuum brightness distribu- 
tion. The central active star-forming regions are detected in all the 
galaxies. The peak fluxes with 1 a errors are listed in Table [T] All 
the sources are more extended than the beam (~ 5 arcsec), and the 
total fluxes are also listed in Table [T] The high brightness region 
(> 2a) is concentrated, considering that the interferometry is sensi- 
tive to an extension of 19 arcsec. Hereafter, we use the term 'central 
star forming-region' to indicate the region where the SMA bright- 
ness is > 2(7. If we compare the SCUBA im ages of II Zw 40 and H e 
2-10 at 450 \xm (beam size ~ 8.5 arcsec) in Galliano et al.l l l2005h . 
the size of the centrally concentrated bright regions are consistent, 
but we miss the diffuse component (shown by the lowest contours 
extending > 1 arcmin in their figure 1). The diffuse component is 
also clear in their 850 [im image with a larger beam (~ 15 arcsec). 
We mainly analyze and discuss the SMA continuum data, although 
we additionally use the CO(3-2) data in Section [T6l 

In Fig. [T] the SMA 880 \im brightness distribution is over- 
laid with the Hubble Space Telescope (HST) Advanced Camera for 
Surveys (ACS) data in the optical (at F550M, F550W, and F550W 
bands for He 2-10, NGC 5253, and II Zw 40, respectively]! The 
positional accuracy of our SMA image is < 0.1 arcsec, while that 
of the HST is .5 arcsec, limited b y the positional uncertainties 
of the guide stars ( lLaskeretal.lll990h . The centre of submillimetre 
emission coincides with that of optical for II Zw 40, while there is 
a significant shift in He 2-10. As shown bv lKobulnickv & JohnsonI 
fl999) , the 3.6 cm radio image also shows a poor correlation with 
the HST optical image for He 2-10. The poor positional coinci- 



^ In this paper, we call the luminosity of the entire galaxy 'global luminos- 
ity'. 
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Table 1. Data. 



Object 


R.A. 


Dec. 






peak flux'^ 


total flux'' 


free-free'^ 


fraction'^ 


12 + log(0/H)f 




(J2000) 


(J2000) 


(Mpc) 


km s^^ 


(mjy/beam) 


(mJy) 


(mJy) 






He 2-10 


081^ 36™ 15! 1 


-26°24'33'.'8 


10.5 


850 


12.2 ±1.6 


26.2 ±2.3 


12.8-15.4 


0.45-0.64 


8.93 


NGC 5253 




-31°38'24'.'4 


3.7 


401 


37.9 ±2.6 


67.0 ±3.3 


48 ±4 


0.63-0.82 


8.14 


II Zw 40 


05'^55™42!6 


+03°23'31'.'8 


10.5 


767 


8.43 ± 1.15 


13.6 ± 2.0 


10-13.5 


0.64-1 


8.13 



Distance. References: CH?) for He 2-10;'Caldwell & Phillips ( 1989) for NGC 5253; HI 1 for II Z w 40. 

^ Velocity relative to the local standard of rest. References: Vanzi et al. (2009) for He 2-10 and NGC 5253: [Sage et alj fT993) for II Zw 40. 
Peak and total continuum fluxes of our SMA 880 i^m observations. The errors are estimated from the 1 cr background fluctuation. 
The contribution from free-free emission at 880 |xm. 

The fraction of free-free contribution^ 

f Oxygen abundance in the literature: iKobulnickv. Kemiicutt. & Pizagnd i 19991) for He 2-10 and NGC 5253: lThuan & Izotovl i2005h for II Zw 40. 



dence can be attributed to the extinction in the optical image. The 
optical image of NGC 5253 has a complicated morphology, and the 
peak of the SMA image does not match the brightest region in the 
optical. The different morphologies between optical and submm (or 
radio) brightnesses, especially for He 2-10 and NGC 5253, indicate 
that the optical light may be tracing completely different regions 
from those observed in submm and radio. 

Submm continuum emission in galaxies is usually dominated 
by dust thermal radiation an d contaminated by free-free emis- 
sion (e.g. lGalliano et alj2q05h. Both thes e emission processes trace 
star-forming regions (e.g. Icondorj|l992l) . In order to separate the 
free-free contribution at 880 \im, we compare our submm fluxes 
with radio interferometric fluxes in the literature as explained in 
the following subsections. We neglect the non-thermal component 
throughout this paper since the radio spectrum at wavelengths < 
a few cm is consistent with free-free emission in the central star- 
forming regions (see each subsection below). 



2.1 He 2-10 

IKobulnickv & JohnsonI f 19991) observed He 2-10 by the Very Large 
Array (VLA). Their 3.6 cm continuum image (with a beam size of 
0.70 arcsec x 0.57 arcsec FWHM) fairly traces our SMA bright- 
ness, although individual knots detected in the VLA image are not 
resolved by our SMA observation. They also have high-resolution 
images at 2 and 6 cm, where the measured fluxes are 17.5±0.7 mJy 
and 14.5±0.3 mJy, respectively. The 3.6 cm flux whose {u, v) cov- 
erage is matched to these maps is 15.6±0.4 mJy. These three fluxes 
are consistent with the flat spectrum expecte d for free-free emis- 
sion (see also IKobulnickv & Johnsod [l999h . Since the spectrum 
at 2-6 cm is flat, a large free-free absorption at this wavelength 
range can be rejected. Therefore, we use the interferometric 2 cm 
flux (17.5 ± 0.7 mJy) to extrapolate the contribution from free-free 
emission at 880 [im. The flux estimated in this way may be a lower 
limit, since these matched data are sensitive to the compact knots 
and the structures surrounding them, but possibly miss the diffuse 
south-east extension in our SMA image (Fig. [T). This extension is 
ju st along the possible tida l tail pointed out by a CO observation 
in IKobulnickv et al] Il995[) . Thus, we also use the total VLA flux 
(21.1 mJy at 2 cm) as an upper limit. If we con vert these fluxes t o 
880 \im fluxes by assuming dependence ( IOsterbroc3ll989l) . 

we obtain 12.8 mJy and 15.4 mJy for the lower and upper limits of 
the free-free contribution, respectively. These values indicate that 
45-64 per cent of the SMA 880 |i.m flux is free-free emission. 



2.2 NGC 5253 

iMeier. Turner. & Beckl ( l2002h show that the radio continuum spec- 
tral slope of NGC 5253 is consistent with free-free emission. The 
slight extension to the west of the peak is common between their 
3. 1 mm image and our SMA 880 \im image. We adopt the t otal flux 
in the inner 20 arcsec region (54 ± 5 mJy at 3.1 mm; Me ier et al] 
|2002|) . to which SMA can be sensitive. The flat spectrum at > 2 
cm show that free-free emission is optically thin at 3.1 mm. We 
estimate the 880 [im free-free flux by assuming a frequency depen- 
dence of for the free-free spectrum. Then we obtain 48 ± 4 
mJy for the contribution from free-free emission at 880 (xm. Thus, 
63-82 per cent of the SMA 880 vim flux is free-free emission. 

2.3 II Zw 40 

This galaxy has already been analyzed and reported in Hll. The 
centimetre emission from the central part can be fitted by the free- 
free flat spectrum. Since the VLA 2 cm observation is only sensitive 
to the structures smaller tha n 4 arcsec, we al so adopt the single 
dish flux at 2 cm (18.5 mJv; lBecket"ai]|2002l) as an upper limit. 
Hll also justifies the wavelength dependence of at < 2 cm 

for free-free emission by spectral fitting. If we adopt 14-18.5 mJy 
for 2 cm flux, we obtain 10-13.5 mJy for the contribution from 
free-free emission at 880 \im. Thus, 64-100 per cent of the flux 
detected at 880 [im by SMA is free-free emission from the central 
star-forming region. 



3 MODELS 

In order to interpret the 880 [im emission in our sample, we need 
to model two major emission mechanisms (thermal free-free radi- 
ation and dust emission) from the star-forming regions. In particu- 
lar, the total luminosity of dust emission (called FIR luminosity) is 
used to investigate the radio-FIR luminosity relation later Physical 
quantities governing free-free and dust emissions are introduced. 
One of the basic quantities is the stellar mass (A/,) formed at the 
current episode of star formation, since the ionizing photon lumi- 
nosity which determines the free-free emission and the UV lumi- 
nosity which contributes to the heating of dust are proportional to 
Af , . To evaluate Af , , the free-free fluxes extrapolated from radio 
observations in Sections 12. 1112.31 are used (Section 13. 2t . The ob- 
tained Af« is later used to estimate the UV luminosity in Section 
13.31 In addition, we need to determine the mass and temperature of 
the dust to estimate the FIR luminosity, which is used to examine 
the radio-FIR relation. Through the modeling, we can obtain the 
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Figure 1. Contour: SMA 880 |im continuum brightness of (a) He 2-10, 
(b) NGC 5253, and (c) II Zw 40. Solid contours are 2 cr, 3 cr, ... (1 cr = 
1.15 mjy beam~^), while dotted contours are —2 and —3 a. The beam is 
shown in the lower right comer Grey scale: HST ACS optical images, at 
F550M, F555W, and F550W bands for He 2-10, NGC 5253, and II Zw 40, 
respectively. 



stellar mass formed in the current starburst episode, the dust mass 
(the dust optical depth), the dust temperature, all of which are basic 
quantities to understand the strength of starburst and the extent of 
dust enrichment. We use the theoretical models that we applied to 
II Zw 40 in our previous paper (Hll). Some detailed assumptions 
that do not affect our results are simplified. Below we briefly sum- 
marize the models used in this paper The same parameter values 
as those in HI 1 are adopted unless otherwise stated. 



3.1 Basic setups for the star formation 

In HI 1, we modeled the star formation rate (SFR) through the free- 
fall time-scale under a given gas density, while in this paper, we 
treat the SFR as a free parameter This is because the SFR is more 
directly connected to the observed luminosity than the free-fall time 
and the gas density. We assume that the SFR is constant as a func- 
tion of time. As shown below, since we only consider young (< 3 
Myr) star-forming regions, the dependence of the luminosity on star 
formation history is weak in the sense that the total stellar luminos- 
ity is simply determined by the total stellar mass formed in the 
current star formation episode. We assume a Salpeter initial mass 
function (IMF) with a stellar mass range of 0.1-100 Mq. 



3.2 Thermal free-free emission 

Free-free emission contains the information of the total stellar mass 
formed in the current starburst episode as modeled below. We use 
the free-free emission as estimated in Sections 12.1112.31 The total 
stellar mass is necessary to estimate the UV luminosity based on 
which the FIR luminosity is modeled in Section [331 Thus, we first 
relate the free-free emission to the total stellar mass. We assume 
that free-free emission is optically thin since we only consider high 
frequencies such as > 15 GHz (see Sections |2. 1H2.3I for further 
justification). 

The thermal free-free luminosity is proportio nal to the num - 
ber of ionizing photons emitted per unit time, Mon l lCondonlll992h . 
We relate the SFR with following Se ction 4.1 of Hll (origi- 
nally, equation 1 of'Hirash ita & Hunll2006l) . In this paper, we adopt 
the ionizing photon luminosity as a function of s tellar mass by tak- 
ing the solar metallicity case in ISchaereil (l2002h . The dependence 
of NioYi on the stellar metallicity is within a factor of 2 in the metal- 
licity range concerned in this paper 

Because we mainly consider young (< 3 Myr) galaxies, the 
decline of luminosities by stellar death is small. Thus, the luminos- 
ity basically reflects the total stellar mass formed up to the current 
age of the system. In Fig.|2l we show the evolution of the free-free 
luminosity {Lv) at 15 GHz normalized to the stellar mass formed 
(integration of the SFR over the time), which is denoted as A/,. Al- 
though we only show the results at 15 GHz, at other frequencies 
can be estimated by assuming a frequency dependence of oc v~^'^ 
as long as the emission is dominated by free-free and is optically 
thin (Sections |2.m2.3l l. We assume a constant SFR, but the follow- 
ing results are not sensitive to the time variation of the SFR under a 
given Mt as long as we consider young (< 3 Myr) ages. We obtain 
the following formula: 



L^lSGHz) 
2.45 X 10" WHz-i 



(1) 



where ^^,(15 GHz) = 4TrD^,U{15 GHz) [/^(15 GHz) is the flux 
at 15 GHz (2 cm) adopted in Sections r2. H42.3l and listed in Table|2l. 
This equation is not valid if there is a significant contribution from 
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Figure 2. Time evolution of tlie free-free luminosity at 15 GHz (2 cm) 
(solid line; left axis for the scale) and the UV luminosity (dashed line; right 
axis for the scale). Both luminosities are normalized to the stellar mass 
formed. For age < 3 Myr, the decline of luminosities by stellar death is 
negligible so that both luminosities normalized to the stellar mass are con- 
stant. 



non-thermal synchrotron emission or if the age is much older than 3 
Myr. In our SMA sample, the contribution from non-thermal emis- 
sion should be small since the observational spectral indexes at cen- 
timetre wavelengths are consistent with thermal free-free (Section 

13. 

In summary, using L„(15 GHz) estimated from fi^{15 GHz) 
in Sections [2. lM2.3l (Table|2t. we obtain the stellar mass formed in 
the current starburst episode (i\f,) by equation lO. The obtained 
M* is used in Section l331 to estimate the UV luminosity. 



3.3 Dust emission 

To give a physical interpretation to the dust emission component 
in the observed 880 \im flux, we need to determine the mass and 
temperature of the dust. The dust mass (the dust optical depth) and 
the dust temperature are basic quantities to understand the strength 
of starburst and the extent of dust enrichment. The total FIR lu- 
minosity can also be derived and used to examine the radio-FIR 
relation (Section [3. 5t . We do not directly use optical and UV ob- 
servations, since they may only trace regions with less extinctions 
and may not be relevant to what we are modeling (Section |2j. We 
use the same simple model as Hll for the dust emission associ- 
ated with the central star-forming region in each individual galaxy. 
Here we summarize the model. At submm wavelengths, the 'large' 
grains, which achieve radiative equilibrium with the ambient stel- 
lar r adiation field, are the dominant component in the luminosity 
(e.g. lGamanoirZll2005h . For simplicity, we assume that the dust 
is distributed in a thin shell at a distance i?dust from the centre and 
that the young stars are located at the centre of the shell (as also 
assumed in Galliano "et al.ll20()5l) . The following results are valid 
as long as the starbursts are more concentrated than the dust dis- 
tribution. In our sample, the stars associated with t he current star- 
burst are more compact than the dust distribution (lGoriianlll99i 
I Johnson et al1l200(]| ; IVanzi et al.ll200 g). In reality, we should also 
note that it is difficult to constrain the spatial distribution of stars 
because of large dust extinctions. 



Under the above assumptions, the dust optical depth, rdust, for 
the radiation from stars is estimated as 

3A/dust 



Tdust = 



(2) 



where A/dust is the total dust mass in the shell, a — 0.1 [im is the 
typical grain radius (we adopt the geometrical cross section for the 
absorption c ross section), and s = 3 g cm""^ is the grain material 
density (e.g. lOraine & Leelll984l) . 

We assume that the UV luminosity (denoted as £ob) is equal 
to the bolometric luminosity of the OB stars (stars heavier than 
3 Mq); that is, the total UV luminosity is calculated by summing 
all the contribution from OB stars under the star formation his- 
tory and the initial mass function given in Section 13. II In Fig. |2] 
we show the evolution of the UV luminos i ty cal culated by the 
model (see Hll for details). Since 'Schaered (|2002|) does not pro- 
vide the stellar luminosities fo r solar metallicity, we adopt th e zero- 
metallicity case. According to lRaiter. Schaerer. & Fosburvl ( l2010h . 
the variation of UV luminosity by metallicity under a fixed SFR is 
within 0.2 dex if the Salpeter IMF is adopted. Following the result 
shown in Fig.|2l we obtain 



Cob = 5.8 x 10^ 



-0, 



(3) 



where M* is already given by the observed radio flux through equa- 
tion l[T). We assume that the OB stars are located at the centre, so 
that we can estimate the total FIR luminosity, -Lpir, by 



(4) 



where Tdust is given by equation The dust temperature, Tdust, 
is related to the dust mass, 



LpiR = 1.09 X IQ-^^A'/dustTlust, 

The monochromatic dust emission flux is expressed as 

fdust{l^) = Ki^Afdust-Bi/(2dust)/-D^, 



(5) 



(6) 



where B^, is the Planck function, is the mass absorption coeffi- 
cient of the dust, and Tdust is the d ust temperature. W e assume that 
= 0.7(i//340 GUzf cm^ g"' jjames et alj2003) . 
The flux from dust /dust (z^) at 880 [im (v = 340 GHz) used in 
equation ([6]l is obtained by subtracting the free-free flux from the 
total SMA flux as described in Section |2] The values of /dust(i') 
at 880 [im are listed in Table |2] Considering the errors, the up- 
per/lower value of the free-free contribution and the lower/upper 
value of the SMA 880 [im flux are used to obtain the lower/upper 
values of /dust- Equations l|5j and ^ indicate that we determine 
the FIR luminosity so that it is consistent with the observed 880 
[im flux as explained at the end of this paragraph (see also Fig. 4 in 
Hll). In HI 1, we gave 7?dust as a free parameter, but in this paper, 
we determine the typical radius for the distribution of the dust asso- 
ciated with the central star-forming region, to which the SMA data 
is sensitive. We provide -Rdust by the intensity-weighted radius as 



r2 

-"•dust 



I{x, y)dxdy 



I{x,y)dxdy, (7) 



where r — {x, y) is the projected position in the image 
and f is the intensity-weighted centre of the image. Although 
our 880 \im images are conta minated with free-free emission, 
iBeck. Turner. & GoriianI (1200 ll) show similar spatial distributions 
between mid-infrared dust emission and radio free-free emission 
in He 2-10, implying that it is reasonable to assume that both dust 
and free-free have similar spatial extent at 880 |i.m at least for this 
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galaxy. For NGC 5253. ' Turner & Beckl ( |2004 show that most of 
the free-free flux comes from the central 1.2 arcsec region, which 
implies that dust emission is more extended than free-free emis- 
sion. In such a case, -Rdust underestimates the actual extension of 
dust. The value of -Rdust for each galaxy is listed in Table |2] After 
all, with -Rdust (given in Table[2j, A/* (derived from the free-free 
luminosity through equation [TJ and Cob (estimated by equation 
O, the unknown parameters are Afdust and Tdust (note that Tdust 
in equation|5]is given by equation|2j. These two unknowns are ob- 
tained by solving equations ^ and 



3.4 Derived quantities 

The quantities derived by the models above (-M, , -Mdust , and Tdust ) 
are listed in Table |2] The stellar mass formed by the current star- 
burst ranges from a few xlO® to ~ 10^ ^^0' and the dust mass 
associated with the star-forming regions is ~ lO^'-lO^ Mq. The 
solar-metallicity sample. He 2-10, has the largest dust content: if 
the dust mass is normalized to the stellar mass, its dust-to- stellar 
mass ratio is 4-10 x 10""^ in comparison with the values 2-3 
xlQ-^ for NGC 5253 and < 5 x 10"^ for II Zw 40. This indi- 
cates that He 2-10 is the most dust-enriched system probably be- 
cause of the highest metallicity. As argued in HI 1 (see also Section 
14. 2b , the dust is preexisting or grown in dense molecular clouds. 
Since dust growth is efficient in high-metallicity environments, the 
highest dust content in the He 2-10 centre among the three sample 
BCDs can be interpreted to be the consequence of the most efficient 
dust growth. Dust temperatu res (~ 50-60 K) hig her than those in 
the Milky Way (~ 15-20 K: iDraine & Le3l 1 9841) are obtained for 
all the sample, supporting that intense star formation is occurring 
in compact regions. Direct constraints on Tdust from FIR data are 
crucial in future high-resolution FIR data around the spectral peak, 
since Lmn is the most sensitive to Tiust than to the other param- 
eters (see equation |5}- However, because -Mdust works as an ad- 
justing factor under a given 880 |J.m flux through equation (|6]l, the 
uncertainty in Lfir is not so large as expected from the uncertainty 
in Tdust- 



3.5 Radio-FIR relation 

In Section |2l we have shown that a large fraction of the SMA 
880 \im. flux in the sample BCDs is contributed from free-free 
emission, while the submm emission is usually d ominated by dust 
on global scales of galaxies jGalliano et al .I2OOA . The large contri- 
bution from free-free emission to the submm emission may cause a 
significant impact on the radio-FIR relation in such a way that the 
radio luminosity is relatively enhanced (HU). Now we examine 
this issue by plotting the radio-FIR relation for the central regions 
in our sample BCDs. 

In Fig. [S] we show the relation between the monochromatic 
luminosity at 15 GHz [-L^(15 GHz)] and the FIR luminosity -Lfir 
for the central star-forming regions in the sample BCDs. -Lfir 
has already been obtained in Section |3^ (see also Table |2j, while 
-L^(15 GHz) is obtained by using the flux at 15 GHz in Tabled 
multiplied by A-kD^ (see also Sections l2.1M2.3b . The upper and 
lower bounds of -Lfir correspond to the lower and upper bounds 
of /dust (the dust flux at 880 iJ.m) given in Table [2] (Recall that 
/dust is the total SMA 880 \im flux minus the free-free contribu- 
tion. If the free-free contribution were not subtracted, we would 
overestimate -Lfir by a factor of 2 or more.) For He 2-10, the up- 
per and lower limits of -L^(15 GHz) are used, for NGC 5253, we 
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Figure 3. Radio-FIR relation for the central stai-forming regions (the di- 
amonds connected by the solid lines for He 2-10 and NGC 5253 and the 
diamond with an arrow for II Zw 40), in terms of the global relations for 
BCDs (squares and crosses). The observational data for the g lobal emis- 
sion from BCDs are taken from lHunt et al.l J2005l) (squares) and lKIein et alj 
(1991) (crosses). We also show the solid line with 515 = 2.80 and the 
dotted lines with 915 = 2.80 ± 0.34 (i.e. ±1 a). 



use the measured values, and for II Zw 40, we plot the lower value 
for L^(15 GHz) to put an upper limit for Lfir (Fig.|3j. 

For comparison, we plot the observational data of global emis- 
sion from BCDs (i.e. the total luminos ity from the entire system) 
in Fig . |3] The sample is taken from iHunt. Bianchi. & Maiolinol 
l!200 5l) for the FIR and 15 GHz global luminosities (squares) and 
Klein . Weiland. & Brinks! fl991) for the FIR and 10.7 GHz global 
luminosities (crosses). The 10.7 GHz luminosity is converted to the 
15 GHz luminosity by assuming the frequency dependence to be 
oc i/^"', where the spectral index (a) is given in iKlein et alj (T99ll) 
as a result of fitting to the data at some available radio frequencies 
(even if we assume (a) = —0.1 for all the sample, the change of 
15 GHz luminosity is too slight to alfect our results). As an ob- 
servational estimat e of the FIR luminos ity, we adopt an empirically 
derived formula bv lNagata et al.l ( l2002l) . who estimate the total dust 
luminosity at A ^ 40 \im by using the IRAS 60 and 100 (im fluxesQ 
The data are summarized in Appendix [A] 

W e adopt the radio-to-FIR ratio as usually used (e.g. lCondonI 
Il992h . Here we define qis as 



gi5 



= log 



Tf 



3.75 X 1012 W 



log 



L415 GHz) 
WHz-i 



(8) 



The average of gis for the global luminosities is ^15 = 2.80 with a 
standard deviation of 0.34. We show the lines with constant qis in 
Fig.|3] We observe that the data points for the central star-forming 
regions in our sample are consistent with the range of 515 (i.e. 
2.80 ± 0.34) that explains the global radio-FIR relation of BCDs. 

Naively, it would be expected that low-metallicity galaxies 
have relatively small amount of dust, so that the stellar emis- 
sion may not be efficiently reprocessed into FIR jHirashita & HuntI 



The contribution at 8^0 |j.m is neglected since the contribution from 
these wavelengths to the total FIR is small (at most 30 per cent). The con- 
tribution at A > 120 |xm is included, while it is not included in the FIR 
luminosity defined bv lHelou. Soifer. & Rowan-Robinsoij il985l) . 
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Table 2. Derived or estimated quantities. 



Object 


^dust 




/"dust 


/,(15GHz) 




A^dust 


Tdust 


f"dust 






(arcsec) 


(pc) 


(raJy) 


(raJy) 


(106 Mq) 


(10* Mq) 


(K) 




(10^ Lq) 


He 2-10 


3.31 


168 


9-16 


21.1-17.5 


11-9.4 


4.7-9.2 


48^6 


0.69-1.4 


3.2-4.0 


NGC 5253 


3.61 


64.8 


12-26 


65 ± 6*^ 


4.3 ±0.4 


0.65-1.5 


57-53 


0.64-1.5 


1.2-1.9 


II Zw 40 


3.28 


167 


< 5.6 


14-18.5 


7.5 ±0.8 


< 3.7 


45<= 


< 0.55 


< 1.8 



The dust flux estimated from the total SMA 880 |x m flux minus the f ree-free contribution (Table[T]. 
° Converted from the 3.1 mm flux (54 ± 5 mjy) in iMeier et alj j2002h by assuming frequency dependence i/ 
The value is for the upper limit of LpiR. 



I2OO8I) . On the other hand, free-free emission does not have such a 
dependence on metallicity (or dust abundance). Thus, we would ex- 
pect systematically smaller gis for low-metallicity galaxies. How- 
ever, NGC 5253 [12 + log(0/H) = 8.14], compared with He 
2-10 [12 + log(0/H) = 8.93], does not follow this expectation, 
which means that even a low metallicity environment can repro- 
cess the stellar light into FIR with a similar efficiency to a solar 
metallicity environment. As we can see in Table |2] the dust opti- 
cal depth is comparable between He 2-10 and NGC 5253, although 
the dust mass is smaller in NGC 5253 than in He 2-10. This is be- 
cause the distribution of dust is more compact in NGC 5253 than in 
He 2-10. A compact geometr y of dust distribution tends to predict 
a high dust temperature (e.g. iTakeuchi et al.llioosi) . explaining the 
observed higher dust temperature in NGC 5253 than in He 2-10. 
The radio-FIR relation of II Zw 40 is still uncertain because of the 
uncertainty in Lpir. However, it is still possible that this galaxy 
has a similar qis to the other two galaxies. If its gis is significantly 
smaller than the other two galaxies, we need to consider the reason 
other than the metallicity, since NGC 5253, which has a similar 
metallicity to II Zw 40, has a gis value as large as that of a solar 
metallicity object, He 2-10. 

The radio-FIR relation of star-forming dwarf galaxies is sim- 
ilar to that of normal g alaxies in spite of the difference in metallic- 



conclude that CO(3-2) is not detected in this galaxy at the expected 
frequency. 

For He 2-10, we derive the CO(3-2) line flux (S'co(3-2))> 
the FWHM of the line in units of velocity (AV), and the 
radius of the CO emitting region [R = 0.7(a6)^/^ where 
a and b are resp ectively the FWHM major and minor axes; 
iMeier et alj ilOOM . For NGC 5253 and II Zw 40, we derive up- 
per limits for the CO(3-2) line flux. The total mass (dynamical 
mass) in the central region traced by SMA can be estimated by 
jMacLaren, Richardson. & Wolfendalelll988l : lMeier et al.ll2002l) 



(9) 



The m olecular mass denoted as Afmoi is estimated as jMeier et"ail 
I2OOII) 



M,„oi = 1.23 X 10* 



Xco \ / 115 GHz 



gal 



)( 



D 



1 Mpc 



Jy km s 



1 I ^'■32/10' 



(10) 



l2()o3 



Wu 



itv (lKleinetalJ |l99ll ; iHopkins. Schulte-Ladbeck. & Drozdovskvl 
et alj[ 



20081). althou gh some specific galaxies show de 



^cogai = Z.J X lu cm (K km s^^)^^ is the Galac- 
tic conversion factor Jstrong et al.l 1988b. Xco is the metallicity - 
dependent conversion factor dArimoto. Sofue. & Tsuiimotolll996i) . 



viations IcannonetaD M. Thus, metallicity (or dust content) CO(3-2)/CO(l-0) line ral,io. We adopt ii; 



32/10 



cannot be the domin ant factor that governs the radio-FIR relation. 
ICarmon et al. I j2005l) find a spatial variation of FIR-to-radio ratio 
by an order of mag nitude in a metal poor dwarf galaxy, IC 2574. 
iDumas et alj ( I2OIII) also show different radio-FIR relations be- 
tween spiral arms and interarm regions in M5 1 . These observations 
also support the above statement that metallicity is not the domi- 
nant factor that varies the FIR-to-radio ratio. The robustness of the 
FIR-to-radio ratio among our BCD sample is further discussed and 
interpreted in Section|4T| 



3.6 Molecular gas 

Gas mass provides the normalization in estimating star formation 
efficiency and dust-to-gas ratio. The observed wavelength range 
also covers the CO(3-2) rotational transition line at 345.796 GHz 
(at rest). Although quite a lot of de tections are reporte d for CO 
in He 2-10 (e.g. Baas, Is rael, & Koom neef 1994; Kobulni ckv et al.l 
ll995l ; IVanzietalj2009l) , we use our data which have the advantage 
of covering the the same (w, v) range. In Table [3] we summarize 
the quantities derived by our CO(3-2) data. 

The spectra around the expected wavelength of CO(3-2) are 
shown in Fig.|4] The emission is only detected in He 2-10. For II 
Zw 40, we failed to obtain the data at > 344.94 GHz. Yet, we can 



0.6 for all the sam ple as a representative value for dwarf starbursts 
jMeier et al .l200lh . The molecular gas mass is compared to the dust 
mass and the stellar mass to obtain, respectively, the dust-to-gas 
mass ratio (© = Mdust/M^oi), and e, [e, = A/./(A/n,oi + M.)], 
which under normal conditions (in particular no molecule destruc- 
tion), provides the star-formation efficiency. 

We detected CO(3-2) only in He 2-10, while we obtained 
only upper limits for NGC 5253 and II Zw 40. Our observations 
are sensitive only to the centrally concentrated component. Dif- 
fuse components, if any, may be resolved out. Table [3] indicates 
that (Afmoi, A'/tot) = (9.4 x lO'^ Mq, 1.0 x 10* Mq) for the 
central star-forming region of He 2-10 with A.V = 40.8 km s~^. 
Ivanzi et al. (2009) reported A V — 5 3 km s"'^ for the central 20 
arcsec in He 2-10. iMeier et al.i ( l200lh estimate (Afmoi, Aftot) = 
(1.4 X 10® Mq, 3.3 X 10* Mq) and (3.7 x lO'' Mq, 9.4 x 
10'' Mq) for He 2-10 and NGC 5253, respectively. Their large 
velocities and masses are due to different spatial scales traced. No 
CO(3-2) detection has been reported for II Zw 40, but the molecu- 
lar mass derived from CO(2-l) and CO(l-O) (> 5 x lO'^ Mq) is 
not contradictory with our upper limit. 

In the He 2-10 centre, Aftot — A/moi, and the stellar mass 
only occupies a small fraction (~ 10 per cent) of the total mass. A 
star formation efficiency of ~ 10 per cen t is near the values derived 
for the Galactic giant molecular clouds jLada. Lombardi. & AlvesI 
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Table 3. Moleculai' gas properties in the center of the sample BCDs. 



Object 


•S'cO{3-2) 


Ay 


R 




A/tot 


Afniol 


^CoZ-'^COgal 


e* 


V 




(Jy kms-1) 


(kms-i) 


(arcsec) 


(pc) 


(10" Mq) 


(lO*^ Mq) 






(10-3) 


He 2-10 


387 ±8 


40.8 ±1.4 


6.48 


330 


100 ±8 


94 ±9 


0.97 


0.08-0.11 


0.50-0.98 


NGC 5253 


< 20=^ 










< 3.8 


6.0 


> 0.53 


> 1.7 


II Zw 40 


< la^* 










< 20 


6.2 


> 0.27 


> 0.27 



3 cr upper limits with AV = 40 km s 



120101) . The dust-to-ga s ratio is signi ficantly smaller than the Galac- 
tic value (~ 6 X 10~'^: [SDitzeil 19781) although the metallicity is sim- 
ilar. This is probably because we only trace the high-temperature 
dust component that is directly heated by the current starburst with- 
out being shielded by other dust co mponents. A similar un deresti- 
mate of dust mass is also reported bv lSun & HirashitaldlOlli) . Thus, 
we may miss a large amount of dust if a small amount of dust effi- 
ciently shields the stellar light. 

For the low-metallicity objects, NGC 5253 and 11 Zw 40, the 
upper limits of molecular gas mass can be used to constrain the 
lower limits for e, and T>. Both objects have significantly larger e, 
than He 2-10. This is interpreted in two ways: (i) the star forma- 
tion efficiency is actually high, or (ii) molecular gas is dissociated 
quickly (< 3 Myr), causing an underestimate of the total gas mass. 
These two possibilities are further discussed in Section |43] 



4 DISCUSSION 

4.1 Free-free contribution at 880 \xm 

In Hll, the 880 |i.m luminosity in the II Zw 40 centre is shown 
to be dominated by free-free emission. This is why Hll suggests 
that free-free dominated emission at ~ 880 \im can be used to 
select extremely young starbursts (< 3 Myr). In this paper, this is 
confirmed in the sense that the central star-forming regions in He 
2-10 and NGC 5253 also have significant contributions of free- 
free emission at 880 [J-m (Section |2]l- There is a hint in Table [T] 
that as metallicity becomes higher the free-free fraction becomes 
lower (i.e. the fraction of dust emission becomes higher). This trend 
should be checked with a larger sample in the future. 

Although free-free emission has a significant contribution 
even at 880 ^m [while global 880 \x.m emission ('global' refers 
to the integrated e mission of the en tire galaxy) is generally domi- 
nated by dust (e.g. ' Hunt et alj|2005l) 1. the radio-FlR relation of the 
central regions in our sample BCDs is surprisingly consistent with 
that defined by using global luminosities as shown in Fig. [3] This 
is probably due to the high dust temperatures, which enhance the 
FIR luminosity and compensates the small contribution from dust 
emission at 880 |i.m. Therefore, it is possible that this compensat- 
ing effect is a key to understand the robustness of the radio-FIR 
relation. 

The robustness of the radio-FIR relation may simply reflect 
the fact that both radio and FIR luminosit ies are good indicators 
of star formation activities ( ICondonlll99i ). If we assume the ro- 
bustness of the radio-FIR relation, we can explain the small con- 
tribution of dust emission at 880 \xm in comparison with free- 
free emission as follows. By assumption, the free-free flux at 880 
[im [/ff (i^)] is just proportional to the total FIR flux, which is 
proportional t o T^„st (for dust mass absorption coefficient oc v^; 
lOraine & Lee|[l984 ). On the other hand, if the Rayleigh-Ieans ap- 
proximation is valid, the 880 (im dust flux, /dust {v), is proportional 



to Tdust- Therefore, fs{v)/fduBt(v) oc T^^^^, which indicates that 
the fraction of free-free emission at 880 (xm is larger for a higher 
dust temperature. 

To examine this temperature dependence of the free-free frac- 
tion, we plot fail') / fd-astiy) at 880 |i.m in terms of Tdust for the 
global emission in the BCDs which were used t o plot the radio- 
FIR relation in F ig.[3](i.e. the samples taken from lHunt et al.l2005l : 
lKleinetalJl99ll) . For these BCDs, we estimate /dust at 880 vim as 
described in Appendix lAl The free-free flux at 880 \im (fg) is es- 
timated by converting the free-free flux at 15 GHz used in Fig.[3]to 
that at 880 \im by assuming a frequency dependence of i^^"'^. We 
also plot our SMA sample for the central parts in BCDs. We find 
that there is a positive trend of fg / /dust with Tdust as expected 
above. We also draw a line with /a //dust oc T^^g^ with an arbi- 
trary normalization in Fig. |5] This line roughly bridges our sample 
and the above samples, although the large scatter of the data im- 
plies that our discussion may be too simplified. Supported by the 
correlation between fg //dust and Tdust, we conclude that the large 
contribution of free-free emission in our observation is a natural 
consequence of picking up active star-forming regions which have 
high dust temperatures. 

The large free-free fraction in the central star-forming re- 
gions of our sample indicates that the subtraction of free-free emis- 
sion at submil limetre wavelengths i s a critical step in studying the 
dust emission. lGallianoetai] ( l2005h find that the submm excess in 
some BCDs needs to be fitted with 'very cold grains', whose de- 
rived abundance may be sensitive to the subtraction of free-free 
emission. Submm ex cess is preferen tially seen in low-metallicity 
dwarf galaxies (Galametz et al.l201 ll) . Since we have assumed that 
the SMA 880 vim flux is dominated by the warm dust component 
heated by young stars, we cannot test the existence of this very cold 
component in the central star-forming regions. We need to wait for 
multi-wavelength observations at FIR-submm wavelengths with 
spatial resolutions of a few arcsec to tackle this issue. 



4.2 Dust abundance 

As mentioned in Section [T4] the dust-to-stellar mass ratio is larger 
in He 2-10 than in the other two galaxies, while as shown in Sec- 
tion |3.6| the dust-to-molecular gas mass ratio does not necessarily 
trace the trend of metallicity. The dust-to-gas ratio derived for He 
2-10 is much lower than the Galactic dust-to-gas ratio (6 x 10"''; 
ISpitzedri978h , although He 2-10 is a solar-metallicity object. As 
interpreted in Section 1X6] the dust mass in He 2-10 may be un- 
derestimated because we only trace the dust directly heated by the 
central star clusters; that is, the dust mass estimate from dust emis- 
sion is not sensitive to the dust in regions where the UV radiation 
from the central stars is shielded. 

As argued in Hll, the observed dust should either be pre- 
existing or have gr own by accretion i n the dense star -forming re- 
gions. According to lHirashita & Kuol j201 ll) (see also llnou3l201ll ; 
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Figure 4. Spectra around the frequency where the CO(3-2) Hne is ex- 
pected. Panels (a), (b), and (c) show the spectrum in the central part of 
He 2-10, NGC 5253, and II Zw 40, respectively. The LSR velocity is also 
shown on the upper axis. For the latter two galaxies, we did not detect signif- 
icant CO(3-2) emission, whose expected frequency is marked by the arrow. 
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Figure 5. Relation between the free-free-to-dust flux ratio estimated at 
880 ^m and the dust temperature. He 2-10, II Zw 40, and NGC 5253, 
in terms of the global relations for BC Ds. The observational data for the 
global emission of BCDs are taken from lHunt et al. (2005) and lKIein et alj 
tl991j) for squares and crosses, respectively. We also show the relation with 



/ff//dus 



where a constant C is chosen arbitrarily. 



lAsano et al .l2012h . the dust-growth time-scale for silicate (a similar 
time-scale is obtained for carbonaceous dust) is estimated as 



2.1 X 10^ yr 

Vl05 



riH 



0.1 \im 



50 K 



IZ, 

-1/2 







Vo.sy 



(11) 



where (a^) and (a^) are the averages of a'^ and (a is the 
grain radius) for the grain size distribution (we adopt (a^) / (a^) = 
0.1 \im}, Z is the metallicity (we adopt Z — 1.73, 0.28, and 
0.28 Zq for He 2-10, NGC 5253, and II Zw 40, respectively, by as- 
suming that the solar oxygen abundance is 12-|-log(0/H) = 8.69; 
lLoddersll2003h . nn is the hydrogen number density (we adopt 
riH = 10^ cm~^ for the dense star-forming regions in BCDs; 
Hll), Tgas is the gas teniperatu re (we adopt Tgas ~ 50 K; 
Iwilson. Walker. & Thomlevlll997h and 5* is the sticking efficiency 
of t he relevant metal species onto the dust surface (we adopt S — 
0.3; lLeitch-Devhn & Williams! 1985l ; lGrassi et al.l20ri ). Then, we 
obtain Tgrow ~ 0.12, 0.75 and 0.75 Myr for He 2-10, NGC 5253, 
and II Zw 40, respectively. These time-scales are shorter than or 
comparable to the ages of the star-forming regions (< a few Myr). 
Therefore, dust growth by the accretion of gas-phase metals can be 
an effective mechanism of increasing the dust content in the central 
star-forming regions of those BCDs. 



4.3 Molecular gas properties 



In Section 13.61 we have shown that the evaluated star for- 
mation efficiencies in the central parts of the low-metallicity 
BCDs (NGC 5253 and II Zw 40) are high. Such high star 
formation efficiencies in low-metallicity galaxies are already 
reported w hen CO emission is used for the molecular gas 
tracer (e.g. ISchruba et al.l l2012h . However, no physical mecha- 
nism that make s the star formation in metal-poor gas efficient 
is known. Hira shita & Ferra ra (2002) suggest that star formation 
is rather inefficient in low metallicity (i.e. low dust-to-gas ra- 
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Figure 6. Relation between FIR luminosity and CO(3-2) luminosity for 
the central part of He 2 -10, NGC 5253, an d II Zw 40. We also show the 
linear fitting relation in iMeier et al ] j200lh applicable for dwarf galaxies 
(solid line) and for non-d warf galaxies (dotted line; originally derived from 
iMauersberger et al. 1999). The dashed lines show the typical scatter for the 
dwarf sample in flVteier et al.. ilOOU) . 



tip) environments becaus e of inefficient H2 foiination (see also 
iGnedin & Kravtso\il201 ll v Another inteipretation of the high star 
formation efficiencies is that CO is more easily dissociated in 



sociating photons by dust 


Wolfire, HoUenbach. & McKedUoid: 


iRrumholz, Lerov, & McKee 


201 1 ; Shettv et alj|201 1 ). Lerov et alj 


(2011) sussest that CO dissociation effect becomes dominant be- 



low 12 + log (O/H) = 8.2-8.4, which is consistent with CO defi- 
ciency of NGC 5253 and II Zw 40. In such a case, CO is not a good 
tracer of the total gas mass, and the gas mass estimated from CO is 
an underestimate. The high dust-to-gas ratio in NGC 5253 can be 
an artifact of an underestimate of the molecular gas mass. 

In order to present the deficiency of CO more clearly, we 
show the relation between CO luminosity and FIR luminosity in 
Fig- HI As a reference, we als o show the g l obal ( i.e. entire galaxy) 
relation of a dwarf samp le in iMeier et all ( 1200 lb (originally from 
IMauersberger et al.l 19991) . Compared with this relation (solid line), 
our SMA data point for NGC 5253 is deviated downwards. This 
indicates that CO is depleted in the centre while the FIR dust emis- 
sion still has a significant intensity there. The existence of the cen- 
tral dust emission demonstrates the existence of gas since gas and 
dust are generally well mixed and coexisting; thus, the natural in- 
terpretation of the deficient CO emission is that CO is highly dis- 
sociated in these low-metallicity galaxies. Another low-metallicity 
galaxy, II Zw 40, may also be deviated downward, but since we 
only obtained an upper limit for the FIR luminosity, II Zw 40 may 
be consistent with the global relation. The relatively rich molecular 
gas content in He 2-10 may be explained by a stronger shielding 
of dissociating radiation due to its rich dust content compared with 
the other two BCDs. This trend of CO abundance with metallicity 
is consistent with theoretical predictions by Kmmholz et al ( 201 1). 

We also show in Fig.|6]a line for non-dwarf galaxies (mostly 
spiral galaxies). Such a correlation between CO and FIR l uminosi- 
ties has been well investigated (e.g. lGao 8i Solomonl2004h . The re- 
lation for the dwarf sample (solid line) is systematically deviated. 
The deviation is even larger if we focus on the central part of NGC 
5253. It is probably because we pick up a region with an intense ra- 



diation field (i.e., strong CO dissociation) if we pick up the central 
region. 



4.4 Implication for high-redshift observations 

High-redshift starbursts a re often traced by observ ations of dust 
continuum and/or CO (e.g.' Michatowski et al .12010). BCDs, which 
host compact star-forming regions with various evolutionary stages 
(i.e. various metallicities), provide good 'laboratories' of high-z 
galaxies (or galaxy evolution), although our sample size is still lim- 
ited. For the He 2-10 centre, the gas mass derived from the CO 
observation through a conventional conversion factor matches the 
total dynamical mass for the He 2-10 centre, supporting that CO is 
a good tracer of the total gas mass in this particular BCD. However, 
as discussed in Section [T6l large values of star formation efficiency 
in low-metallicity BCDs, NGC 5253 and II Zw 40, imply that CO 
luminosity is not necessarily a good indicator of the total gas mass 
in low-metallicity starburst environments. Thus, if a high-2: galaxy 
has a metallicity significantly lower than solar, CO luminosity has 
the risk of underestimating the total gas mass. 

This underestimate of gas mass propagates to the estimate of 
dust-to-gas ratio. Moreover, there is a potential risk of underesti- 
mating dust mass because we may miss low-temperature dust com- 
ponents, which are less luminous than high-temperature compo- 
nents directly heated by the current starbursts (Section [4. 2| l. Thus, 
as long as we trace the dust content by emission and the gas content 
by CO lines, there is a possibility that the obtained dust-to-gas ratio 
is not reliable. 

For the origin of dust, it is suggested that dust growth is the 
most e fficient mechani s m of dust productio n in high-redshift star- 
bursts ( lMattssonll20Tll ; IValiante et al.ll2oT ij). As suggested in our 
BCD observations (Section [4.2) , high gas densities in intense star- 
bursts make the grain growth time-scale short. 



5 CONCLUSION 

In order to reveal the submm-radio radiative properties of young 
active starbursts in BCDs, the central star-forming regions in He 
2-10, NGC 5253, and II Zw 40 were observed in the 340 GHz 
(880 \im) band at ~ 5" resolution with SMA. The 880 |J.m fluxes 
have been decomposed into free-free and dust emission compo- 
nents by using centimetre radio data in the literature. At 880 [im, 
free-free emission has proven to have a contribution comparable 
to or larger than dust emission in the central part of all the three 
BCDs. We have also shown that the fraction of free-free emission 
at 880 \im has a positive correlation with the dust temperature. In 
spite of the dominance of free-free emission at 880 vim, the radio- 
FIR relation of the central parts in the BCDs is consistent with the 
relation defined for the global luminosities (i.e. the luminosities 
in the entire system), supporting the robustness of the radio-FIR 
correlation. Finally, we have analyzed the CO(3-2) emission line, 
finding that CO is a good tracer of the total gas mass in He 2-10, 
while CO is deficient in the low-metallicity BCDs (NGC 5253 and 
II Zw 40), probably due to inefficient shielding of dissociating pho- 
tons in dust-poor environments. We also point out that the dust mass 
may be potentially underestimated since emission is always biased 
to the dust directly heated by the stars. Thus, we should keep in 
mind these potential underestimates of gas and dust masses in in- 
terpreting submm observations of not only nearby galaxies but also 
high-redshift galaxies. 
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APPENDIX A: DATA FOR THE RADIO-FIR RELATION 

We collected radio fluxes around 2 cm (15 GHz) and FI R (IRAS 60 
and 10 yi m) fluxes o f BCDs in the literature. We adopt Hunt e t alj 
( l2005h and lKleinet al. (1991) (Tables [ATI and [A2l respectively). 
The flux listed in these two tables is a g lobal flux (i . e. flux from the 
entire galaxy). For the radio fluxes in iHunt et al.l (l2005h . if more 
than two data are available at the same wavelength, we adopt the 
data in the lowest-resolution mode, which is expected to have the 
smallest missing flux. We select BCDs with de tections both at IRA S 
60 and 100 ^lm, and at 2 cm (15 GHz). For lKlein etZI jl99lh 's 
sample, we adopt the Effelsberg 2.8 cm flux and convert it to the 
2 cm flux by adopting average radio spectral indices (a) given in 
Table 6 of lKlein et alTjl99lh. We exclude BCD s without detection 
by I RAS. The /j?A5 fluxes of lKlein et alj Il99lh 's samples are taken 
from lMoshir et al.l ( Il990l) . 

For these samples, we need to estimate dust temperatures and 
fluxes at 880 [xm in Fig. [5] We adopt the follo wing estimator, rLG2, 



for the dust temperature ( Nagata et al.ll2002l) 



T.G. = 11.8 + 13.8 [K]. (Al) 

We use this dust temperature for Tdust in the text. By using Tlg2, 
the dust flux at 880 [im is estimated as 

/duBt(^)l8S0..n - (100/880) ,,p(ig _ 1 

x/^(100nm). (A2) 



This paper has been typeset from a TgX/ MTgX file prepared by the 
author. 
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Table Al. lHunt etai] j2005l) 's sample for the radio-FIR relation. 



Galaxy 


M2 cm) 


/.(60^m) 


/^(lOO ^m) 


D 




(mJy) 


(Jy) 


(Jy) 


(Mpc) 


He 2-10 


21.1 ± 1.2 


24.1 ±2.4 


26.4 ±2.6 


10.5 


NGC 5253 


54 ±5 


30.5 ±1.2 


29.4 ± 1.8 


3.7 


II Zw 40 


12 ±3 


6.61 ±0.70 


5.80 ±0.90 


10.5 


Mrk33 


16.0 ±0.06 


4.77 ±0.04 


5.99 ± 0.13 


24.9=^ 



'^ lTull-vUl988l) . 



lable A2.IKlein 


et alj (19911) s 


sample for the radio- 


-FIR relation. 






Galaxy^ 


/.(2.8 cm) 


(a) 


/,(60nm) 


/.(lOO^m) 


D 




(mJy) 




(Jy) 


(Jy) 


(Mpc) 


Haro 14 


2.4 ±0.8 


-0.51 ±0.18'' 


0.530 + 0.069 


1.04 + 0.14 


13.9 


Haro 15 


6.6 ±0.5 


-0.83 ±0.21 


1.35 + 0.12 


1.97 + 0.20 


95.4 


Mrk 370 


2.7 ±0.7 


-0.44 ±0.11 


1.21 + 0.11 


3.03 + 0.24 


12.9 


II Zw 40 


21 ± 2 


-0.20 ±0.05 


6.61 + 0.70 


5.80 + 0.90 


10.5 


Haro 1 


21 ±3 


-0.48 ±0.08 


8.57 + 0.43 


12.9 + 0.65 


52.0 


Mrk 86 


7.4 ±2.0 


-0.30 ±0.45 


3.24 + 0.19 


6.45 + 0.39 


7.0 


Mrk 401 


11±4 


-0.51 ± 0.18 


2.57 + 0.13 


4.01 + 0.24 


22.2 


Haro 23 


5.5 ±0.8 


+0.71 ±0.30 


0.401 + 0.044 


0.778 + 0.132 


17.4 


Mrk 140 


18 ±6 


-0.29 ±0.04 


0.370 + 0.041 


0.630 + 0.126 


22.7 


Haro 2 


7.1 ±0.9 


-0.59 ±0.10 


4.68 + 0.28 


5.32 + 0.32 


20.5 


Haro 3 


9.1 ±2.0 


-0.25 ±0.02 


4.95 + 0.40 


6.75 + 0.41 


13.9 


Mrk 186 


7±3 


+0.44 ±0.54 


1.09 + 0.066 


2.52 + 0.123 


11. 1 


Mrk 169 


8±2 


-0.49 ±0.07 


3.17 + 0.29 


4.79 + 0.29 


17.4 


Haro 28 


1.5 ±0.7 


-1.11 + 0.60 


1.09 + 0.08 


2.29 + 0.23 


10.7 


Mrk 49 


4.5 ± 1.5 


-0.46 ±0.52 


0.724 + 0.065 


0.823 + 0.140 


18.2 


Haro 9 


17 ±5 


+0.33 + 0.41 


2.63 + 0.21 


4.47 + 0.27 


13.9 


Mrk 59 


6.1 ± 1.3 


-0.43 + 0.09 


1.97 + 0.12 


2.46 + 0.20 


11.6 


II Zw 70 


4.4 ± 1.0 


-0.26 + 0.06 


0.714 + 0.050 


1.24 + 0.12 


17.1 


Mrk 297 


22 ±5 


-0.81 + 0.09 


6.15 + 0.31 


10.2 + 0.5 


63.0 


Mrk 313 


7.5 ±0.5 


-0.63 + 0.13 


3.80 + 0.30 


7.40 + 0.59 


30.8 


Mrk 314 


3.9 ±0.7== 




1.25 + 0.10 


1.49 + 0.37 


31.1 


III Zw 102 


16.7 ± 1.0 


-0.62 + 0.01 


9.33 + 0.56 


17.8 + 1.1 


25.0 



II Zw 40 and Haro 2 (= Mrk 33) are also included in lHunt et al] j2005l)'s sample (Ta blelATl. 
^ Since (a) is not available, we assume t h e mea n value obtained in lKlein et alj jl99l]) . 

Since the flux given by iKlein et alj Jl99lh is probably overestimated because of confusion with 
nearby sources, we replace d the value with an interferometric flux at 2 cm (15 GHz) given by 
'Peeg. D uric. & Brinks! fi997l) . 

Note: Mrk 527 was originally included in the sample of lKlein et alj il99lh . Because of the suspected con- 
tamination of nearby sources, this galaxy is not included in our sample. 
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